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ABSTRACT 

O I Both diffuse high energy gamma-rays and an extended electron-positron an- 

^ I nihilation hne emission have been observed in the Galactic Center (GC) region. 

^ I Although X-ray observations indicate that the galactic black hole Sgr A* is in- 

active now, we suggest that Sgr A* can become active when a captured star is 
^ _ tidally disrupted and matter is accreted into the black hole. As a consequence 

^ I the galactic black hole could be a powerful source of relativistic protons. We 

are able to explain the current observed diffuse gamma-rays and the very de- 
>■ ' tailed 511 keV annihilation line of secondary positrons hj p — p collisions of such 

protons, with appropriate injection times and energy. Relativistic protons could 



in 



^ ' have been injected into the ambient material if the black hole captured a SOM© 

■ star at several tens million years ago. An alternative possibility is that the black 



O 



hole continues to capture stars with ~1M0 every hundred thousand years. Sec- 
ondary positrons produced hj p — p collisions at energies > 30 MeV are cooled 
Qh' down to thermal energies by Coulomb collisions, and annihilate in the warm 

O ' neutral and ionized phases of the interstellar medium with temperatures about 

+-* ' several eV, because the annihilation cross-section reaches its maximum at these 

■ temperatures. It takes about ten million years for the positrons to cool down 
to thermal temperatures so they can diffuse into a very large extended region 

^ ' around the Galactic center. A much more recent star capture may be also able 

■ to account for recent TeV observations within 10 pc of the galactic center as well 
as for the unidentified GeV gamma-ray sources found by EGRET at GC. The 
spectral difference between the GeV flux and the TeV flux could be explained 
naturally in this model as well. 
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- black hole - radiation mechanisms : nonthermal 
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1. Introduction 

The annihilation of cosmic ray positrons produced as a result of cosmic ray proton 
collisions with the ambient plasma was discussed in a number of papers (see, e.g., Ginzburg 
and Syrovatskii, 1964; Hayakawa et al. 1964; Stecker 1969; Bussard et al. 1979, etc.). A 
typical cosmic ray positron resulting from secondary 7r-meson decay may undergo one of three 
fates: 1) escape from the Galaxy, 2) annihilate with an electron while at relativistic energy 
{in-flight annihilation), or 3) lose almost all its energy before annihilation. One can consider 
that annihilation may occur either between free electrons and positrons {free annihilation), 
or through the formation of the intermediate bound state of positronium. Positronium can 
form in either a singlet state, which annihilate into two 511 keV photons, or in a triplet state 
which decays by three photon annihilation producing a continuum emission below 511 keV 
(see Ore and Powell 1949). 

The cross-section for positron annihilation as a function of energy was given first by 
Dirac (1930) (see also Heitler 1960). 

Radiation from electron-positron annihilation from the central region of our Galaxy was 
first reported more than 30 years ago (Johnson et al. 1972). In 1978 it was reported that there 
is a flux of the positron annihilation line (see Leventhal et al. 1978) which might be produced 
in a compact source at or near the Galactic center. First estimations of the morphology 
of this fine was obtained by the OSSE telescope on-board the Gompton observatory (see, 
e.g., Gheng et al. 1997). With the launch of the INTEGRAL telescope high resolution 
spectroscopy with reasonable angular resolutions became available. These measurements 
showed a flux of the annihilation emission from the central bulge. The bulge annihilation 
emission is highly symmetric with an extension of 8°. The INTEGRAL team measured 
a 511 keV line, of width 2.2 keV, with a flux (1.01 ± 0.02) x 10"^ ph cm^^g-i ^nd an 
ortho-positronium annihilation continuum flux equal TO (4.3 ± 0.3) x 10^'^ ph cm^^s^^ (see 
Teegarden et al. 2005; Knoedlseder et al. 2005; Ghurazov et al. 2005; Jean et al. 2005). 

Prom the ratio between line and continuum intensities near and below 511 keV it was 
concluded that positrons annihilate via positronium formation. The fraction of positronium 
was estimated to be in the region of 92-97% (Weidenspointner et al. 2006; Jean et al. 2005). 
From the observed width of the annihilation line it was concluded (see Ghurazov et al. 2005; 
Jean et al. 2005) that the annihilation of thermal positrons takes place in a relatively warm 
(several eV) and lowly ionized (~ 0.1) interstellar medium. 

The origin of such positrons is still poorly understood. A large variety of positron sources 
have been proposed. Among them novae and supernovae stars (Ghan and Lingenfelter 1993; 
Dermer and Murphy, 2001). However, as follows from Weidenspointner et al. (2005) the 
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annihilation of positrons appears to be even more concentrated in the bulge than are old 
stellar populations such as Type la supernovae, novae, or low-mass X-ray binaries. New, 
speculative, physics such as positron production from hypernovae stars (Casse et al. 2004), 
light dark matter (see Boehm et al. 2004) etc. has begun to be discussed as a possible 
solution. A supermassive black hole of 2.45 x 10® Mq is another candidate as a source of 
positrons. Thus, Zheleznyakov and Belyanin (1994) concluded that the annihilation emission 
is consistent with a model in which electron-positron pairs, created in the vicinity of a black 
hole, are trapped in the magnetosphere supported by an accretion disk. 

Secondary positrons generated as a result of cosmic ray collisions with the ambient 
plasma is another alternative model. Mclia et al. (1998); Fatuzzo et al. (2001) and Fatuzzo 
and Melia (2003) suggested a supernova coinciding with the central radio source Sgr A East 
as an emitter of protons which can produce gamma-rays as well as high energy secondary 
positrons. However, they showed from simplified equations that Sgr A East could not be the 
source of annihilation radiation from the Galactic center because, based on their estimates, 
the thermalization time was too long. Besides, the annihilation emission appears to be 
diffuse. So far, there is no evidence for emission from point-like sources (Knoedlseder et 
al. 2005; Bouchet et al. 2005). On the other hand, the EGRET telescope found a flux of 
gamma-rays from the Galactic center at energies > 500 MeV (Mayer-Hasselwander et al. 
1998) which was seen as a strong excess of gamma-rays peaking in an error circle of 0.2" 
radius surrounded by a strong emission maximum within ~ 5 — 8° degrees from the Galactic 
center. 

In our model we assume the black hole is a source of high energy protons generated by 
star accretion. They produce secondary gamma-rays and relativistic positrons as a result 
oi p — p collisions. The secondary positrons requires over 10 million years to cool down 
before annihilation. Hence they are able to propagate far away from the central source 
during their lifetime and to fill a sphere with the radius about several hundred pc. Therefore 
the annihilation emission can be seen as diffuse if these protons were ejected a relatively 
long time ago. From kinetic equations we shall show that processes of Coulomb collisions 
are effective enough to cool down these relativistic positrons and to thermalize them before 
their annihilation, which can explain the origin of the annihilation emission from the Galactic 
center. 



2. Energy Release Supplied by a Black Hole 

We will first summarize how much energy must be released from the Galactic black hole 
in order to explain the observed X-rays, diffuse gamma-rays, the annihilation emission and 
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TeV emission. 



2.1. 



X-rays 



The ASCA telescope found hard X-ray emission from in the energy range 2-10 keV the 
central part of the Galactic plane (Koyama et al. 1996). The total thermal energy of the 
hot plasma component inside a 1° region around the Galactic center was estimated to be 
10^^ — 10^^ erg. A large energy generation rate in the Galactic Center is required ~ 10^^ — 10^^ 
erg s~^ and a supernova origin of this emission seems to be implausible. 

Recently, Chrandra observations showed intense X-ray emission at energy E'-^ ~ 8 keV 
from the inner 20 pc of the Galaxy with a luminosity ~ 10^^ erg s~^. We know of no class 
of object which can provide this luminosity there. Therefore this emission was interpreted as 
thermal emission of an optically thin plasma with a temperature of 8 keV. However, for the 
gravitational potential near the Galactic center (see Breitschwerdt et al. 1993) the escape 
velocity is about 900 km s^^. The sound speed of an 8 keV plasma is 1500 km s^^ and 
the required power necessary to compensate the plasma outflow is ~ 10^° erg s^^. This is 
equivalent to the entire of supernova kinetic energy output occurring every 3,000 years, and 
is unreasonably high (Muno et al. 2004). 



The EGRET telescope (sec Maycr-Hassclwandcr et al. (1998)) found a gamma-ray flux 
toward the Galactic Center of the order of 2 x 10'^'' erg s~^ for energies > 500 MeV in 
an error circle of 0.2" radius. The photon spectrum from the center (ph cm~^s~^) can be 
represented as 



2.2. 



Gamma-rays. 




(1) 



for < 1900 MeV 



and 





for E^ > 1900 MeV 
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The bulk of the emission is most compatible with it originating in a volume with a 
radius 85 pc though, in view of uncertainties, the source excess is marginally compatible 
with emission from a single compact object. If we assume that the central gamma-ray 
source is due to injection of relativistic protons we can estimate their energy. The flux of 
gamma-rays from the central source is about ~ 2 - 10^^ erg s~^. The lifetime of relativistic 
protons in the interstellar medium with a density n ^ 1 — 10 cm~'^ is about Tp ~ 10^"^ — lO^^s, 
which gives a total energy in relativistic protons Wp c:^ ■ Tp 10^^ erg. This is more than 
can be supplied by the most powerful known galactic source of energy, such as a supernova, 
whose total energy release is about 10^^ erg. 

2.3. Annihilation line emission 

On the basis of the reported flux of annihilation line emission (see Teegarden et al. 2005; 
Knoedlseder et al. 2005; Churazov et al. 2005; Jean et al. 2005), the annihilation rate of 
positrons is roughly ~ lO^^s"^. Assuming that each interaction of a relativistic proton with 
an average energy about 1 GeV produces one secondary positron, we find that the injected 
energy rate of protons is Wp ~ 10^° erg s~^. The total amount of injected proton energy must 
be of the order (or smaller than) Wp ~ WpTc ~ lO^^erg, where Tc ~ 10''^^s is the cooling time 
of relativistic positrons in the interstellar medium with a gas density n ~ 1 cm~^. Such a 
huge amount of energy is the most important constraint of the model. In section 2.5, we will 
discuss some possible energy injection mechanisms. 

2.4. TeV gamma-rays 

TeV gamma-ray emission from the Galactic black hole Sgr A* was first reported by 
the Whipple group (Buckley et al. 1997). More recently it has been confirmed by three 
independent groups, Whipple (Kosack et al. 2004), CANGAROO (Tsuchiya et al. 2004), 
and HESS (Aharonian et al. 2004) at a luminosity > 10^^ ergs s~^. If these TeV gamma- 
rays are produced via p-p collisions, the possible input proton energy could be 10^°~^^er(7 
depending on the proton injection spectrum and the diffusion coefficient (Aharonian and 
Neronov 2005a). 

Thus, observational data indicate a huge energy release at the Galactic center with an 
energy of the order > 10^^ erg, which cannot be supphed by known Galactic sources, and a 
significant part of this energy is enclosed in the form of relativistic protons. One can assume 
that such energy can be released when a star is disrupted near a black hole. 
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2.5. 



Theoretical estimates. 



The rate at which a massive black hole in a dense star cluster tidally disrupts and 
swallows stars has been studied extensively (e.g. Hills 1975; Bahcall and Wolf 1976; Lightman 
and Shapiro 1977). Basically when a star trajectory happens to be sufficiently close to a 
massive black hole, the star would be captured and eventually disrupted by tidal forces. After 
a dynamical time-scale (orbital time-scale) , the debris of a tidally disrupted star will form a 
transient accretion disk around the massive black hole, with a radius typically comparable to 
the tidal capture radius (Rees 1988). Rees has also argued that most of the debris material 
will be swallowed by a black hole with a mass ~ IO^Mq on a time scale of ~ 1 yr for a 
thick hot ring, or ~ 10^ yrs for a thin cool disk, respectively. The capture rate is essentially 
a problem of loss-cone diffusion-diffusion in angular momentum rather than energy. By 
assuming a Salpeter mass function for the stars, Syer and Ulmer (1999) have estimated the 
capture rate in our Galaxy as ~ 4.8 x 10~^yr~^ for main sequence stars and ~ 8.5 x lQ~^yr~^ 
for red giant stars, respectively. However, the actual capture rate depends sensitively on the 
assumed mass function of stars, the stellar evolution model used, the radius and mass of the 
captured star, the black hole mass and the internal dispersion velocity of stars {vg) around 
the black hole. For example, Cheng and Lu (2001) obtained a longer capture time ~ 10^ 
years, by taking Vg = lO^km/s and M^h = 2.45 x IO^Mq. Therefore the capture time for a 
main sequence star with mass ~ IM© could range from several ten thousand years to several 
hundred thousand years. The capture time for the more massive stars is expected to be even 
longer. For t> tpeak, the accretion rate evolves as, Rees (1988); Phinney (1989) 



where and are the mass and the radius of the captured star, respectively and tpeak ~ 
1.59tmin, where 



is the characteristic time for the debris to return to the pericenter (Lu et al. 2005). The recent 
Chandra observations of three large amphtude, high-luminosity soft X-ray flares in AGNs 
provide strong evidence for the tidal capture events, and the decrease of X-ray luminosity 
indeed follows closely the above theoretical predictions (e.g. Halpern, Gezari and Komossa 



Sirota et al. (2005) analyzed the tidal disruption of stars near a black hole, a process 
which is a possible mechanism for the activity of galactic nuclei. They concluded that the 
total power released during the destructions of a star is on average 4 • 10^^ erg s~^. However, 
the duration of the energy emission due to the absorption of a star (~ 200 years) is much 




(3) 





2004). 
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shorter than the mean time between two absorption events. Therefore most of the time 
galactic nuclei are in a low active state. Furthermore, the maximum energy release must be 
much higher than the mean estimate. For example, the highest luminosity observed in PG 
0052+251 (see Table 1 in Sirota et al. 2005) is 5 x 10^"^ erg s~^, which gives a total energy 
release during 200 years of about 3 x lO^'' erg. 

It is very important to know how much of the accretion power will be converted into 
the out-fiow of relativistic protons. Processes of particle acceleration near black holes are 
not well known though several models of particle acceleration in accretion disks and jets of 
black holes have been developed (sec, e.g., Kardashev 2001; Heinz and Sunyaev 2002; Le 
Truong and Becker 2004; Aharonian and Neronov 2005b). For our aims we estimate roughly 
from Eq.(3) the energy of disruption which can be transferred to protons. In the case of star 
capture, it is very natural to assume that the resulting jet contains mainly protons simply 
because of the dominance of hydrogen in stars. We do not know how the jet is formed but if 
we take the analogy of Active Galactic Nuclei (AGN), then we can argue that the observed 
Doppler factor of AGN is of order 10. and we just assume that the jet coming out from our 
galactic black hole is mildly relativistic and its main composition is protons. 

Falcke and Biermann (1999) have argued that the conversion efficiency (rjp) from accre- 
tion power (Mc^) into the the energy of jet motion ranges from 10"''" to 10~^. Integrating 
the Eq.(3) and using typical values of the parameters, the energy carried away by relativistic 
protons is estimated as 

A^;^ ~ 6 X 10^2(VlO-')(M,/MQ)erg. (5) 

In addition to the accretion power, Cheng and Lu (2001) have argued that if the transient 
accretion disk can generate a sufficiently strong magnetic field, due to the instability of the 
disk, this strong magnetic field can initiate the Blandford-Znajek process (Blandford and 
Znajek 1977) to extract rotation energy from the black hole. They estimate that the 
maximum energy that can be extracted from a black hole is given by 

AE^,, ~ 3 X 10'^A^f{A)MieTg. (6) 

where A is the dimensionless angular momentum of the black hole, f{A) is a constant for 
given A and Mg is the black hole mass in units of 10^ Mq. If the black hole is rotating in 
maximum angular velocity, A = 1, then /(74)=1.14. 

The maximum energy in relativistic protons can be estimated from Eq.(5) or (6). If a 
star with the mass about 50 Mq is captured by a black hole, it gives an energy in relativistic 
protons as high as ~ 10^"^ erg. 
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Thus, we conclude that when eventually a massive star is captured, a huge amount of 
energy can be released in the form of relativistic protons during a very short time. Below we 
assume that such primary protons interact with the medium gas in the region, its intensity 
derived from the gamma-ray data, and produce there secondary gamma-rays and positrons. 

In the summary of these two sections we present the conditions which our model of 
annihilation emission by the black hole should satisfy: 

1. The rate of production of thermal positrons at present should be ~ 10^^ s~^; 

2. If positrons are produced via p — p collisions the total energy of relativistic protons 
must not exceed 10^^ — 10^^ erg; 

3. The eruption time of protons is much smaller than any other characteristic time of 
the problem. Therefore the production function of protons can be described as a delta 
function in time; 

4. Though positrons are emitted by the black hole - a point-like source- their spatial dis- 
tribution when they cool down to thermal energies should extend over a region with 
an angular radius 5" — 8° around the Galactic Center; 

5. If the gamma-ray emission near the Galactic center and the annihilation emission have 
a common origin the flux of gamma-rays with energies > 500 MeV should be about 

~ 2 • 10^^ erg s^^ 



3. Proton Spectrum 

We take the source function of protons as 

Q{r,Ep,t)^A{E,)S{r)S{t) (7) 



Particles can be accelerated near the hole by a shock wave appearing in the accreting 
matter. According to typical shock acceleration models the momentum spectrum of accel- 
erated particles is a power- law, with a spectral index 70 between 2-3. If we transform the 
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momentum spectrum to an energy spectrum, then it has the form 

In fitting the observed diffuse gamma-ray spectrum we should be able to provide some 
constraints on 70. However, the error bars of the observed data are quite large (cf. Fig. 
3), and the possible range of 70 is around 2.6-3.1. In fact the best fit gamma-ray spectral 
index is 3.1 for photon energies larger than 2 GeV (cf. Eq.2). However, the exact value of 
7o affects neither the gamma-ray luminosity nor the annihilation flux signiflcantly, because 
most of the energy comes from protons with energies of several GeV. Below we take 70 ~ 2.6 
as reference and the constant Aq will be found from the gamma-ray data. 

The spatial distribution of the protons can easily be derived from the well-known equa- 
tion of cosmic ray propagation (see Berezinskn et al. 1990), 

?^-ViD,VnM§i{fn,)^"^^Q(r.E,,t) (9) 

Here dE/dt is the term describing particle energy losses and Tp is the characteristic time of 
catastrophic losses. 

In the relativistic energy range the lifetime of protons is determined hy p — p collisions, 
and equals Tp ~ 10^^ s for a gas density n = 1 cm^'^ . The rate of continuous energy loss is 
negligible in this energy range {dE/dt ~ 0). The solution of Eq. (9) in this case is similar 
to the solution of the equation of thermal conductivity 



exp 



np{r, Ep, t) = A exp [-t/ Tp] ^^^^ ^^^^^ (10) 
and describes an almost uniform proton distribution within the sphere of the radius r < 

In the non-relativistic energy range the injection spectrum is transformed by ionization 
losses which for a medium density n = 1 cm^^ are significant at relatively low proton energies. 
The rate of ionization losses is (Hayakawa 1964; Ginzburg 1989) 

dE 27re^n / mVM^\ a_ ^^^^ 



dt mc(3{E) \ Ane'^h^n J -/b 

where Wmax is the highest energy transmitted to an ambient electron, and P{E) — v/c. The 
characteristic time of ionization (Coulomb) losses is Tc ~ -B/ (dE/dt). Variations of ionization 
loss times for electrons (dashed-dotted line) and protons (sohd line) are shown in Fig.l. 
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The equation for the total number of particles N{E,t) = Att J np{E,r,t)r'^dr can be 



obtained by integration of Eq.(9) over the volume occupied by the particles when dE/dt ^ 
and is described by Eq.(ll), 



Np{Ep,t)^-- 



'-f+Erf' 



M + £;3/2y^o+2)/3 



1 (70+l)/2 



exp 



^ (12) 



The total spectrum of protons for an arbitrary constant Aq is shown in Fig. 2. Wc sec that 
at the initial stage the spectrum is evolved under the influence of ionization loss but after 
the time t > ~ 10^^ s the p — p nuclear reactions become also essential. 



4. The Gas Components and Particle Propagation in the Bulge Region 

The gas content in the Galactic center is not well known though one can find its char- 
acteristics in Jean et al. (2005). For the interstellar medium several components of the gas 
were suggested from radio and X-ray data (see the review of Ferriere 2001). They are: a cold 
neutral component {n = 10 cm"^, T = 80 K), a warm neutral component {n = 0.3 cm^^, 
T = 8000 K), a warm ionized component (n = 0.17 cm"^, T = 8000 K) and a hot component 
(n = 3 ■ 10~3 cm~3, T = 4.5 • 10^ K). The distribution of these components in the bulge 
region is strongly nonuniform which complicates matters. The point is that gamma-rays 
from TT^ can be produced by all components of the gas if high energy protons penetrate there 
while annihilation photons are produced mainly in the neutral and ionized warm compo- 
nents where the cross-section for annihilation increases by several orders of magnitude as 
compared with the other components. We do not know whether the propagation of high 
energy protons is the same in all components of the bulge medium. We also do not know 
the characteristics of secondary positron propagation: whether they are able to penetrate 
into dense molecular clouds or what is their diffusion coefficient. 

The Bulge (i.e. the region inside a radius ~ 230 pc with a height 45 pc) contains 
7 ■ 10^ Mq of hydro gen gas. 90% of this mass is trapped in small high density (~ 10^ cm ^) 
clouds that gives the average gas density in the bulge central region as high as ~ 700 cm"^ 
while the remaining ~10% is homogeneously distributed with the average density 10 cm^^. 
The rest of the gas in the bulge is contained in a region (like an ellipsoid) with a radius 1.75 
kpc. The HI mass is equally distribute between cold and warm neutral gas. The mass of 
warm and hot ionized gas is about 4-10^ Mq. 90% of this mass is in the warm phase and 
10% is in the hot phase. 
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The diffusion coefficient is unknown in tfie Bulge region. Berezinskii et al. (1990) esti- 
mated tlie average value of the diffusion coefficient in the Galactic Disk from analyses of the 
cosmic ray flux and obtained a value ~ 10^^ cm^s~^. Jean et al. (2005) presented a very 
rough estimate of the diffusion coefficient in the Bulge region using the equation derived by 
Melrose (1980) 

Dc^dJ^^X (13) 

where Db = ^^lv is the Bohm diffusion coefficient, tl and v are the particle gyroradius 
and velocity, respectively, X^ax — 100 pc is the maximum scale of turbulence, 5 — 5/3 for 
a Kolmogorov turbulence spectrum, and 77 ~ ~ 1 is the pressure ratio of magnetic 

fluctuations to the main field. 

For a magnetic field strength B ~ 10~^G (see, e.g. LaRosa et al. 2005) the diffusion 
coefficient of 30 MeV positron is about lO^^cm^s"^. Though the diffusion coefficient decreases 
with particle energy, just this value determines positron propagation in the Bulge region 
because the characteristic time Tc oc E^, i.e. positrons spend most of their lifetime with high 
energy. 

The problem of particle penetration into molecular clouds was analyzed by Skilling and 
Strong (1976); Dogiel and Sharov (1985). In the absence of waves the cosmic ray (CR) flux 
into clouds should be 

j{E) ~ N{E)v (14) 

where v is the particle velocity. However, CRs of a few MeV excite waves to such a high 
level that the flux reduces to 

j{E) ~ N{E)Va (15) 

where the Alfvenic velocity Va « v. 

Therefore, Alfven waves generated outside molecular clouds by the flux of cosmic rays 
are effective at excluding cosmic rays below a few hundred MeV, while particles with higher 
energies freely penetrate into the clouds. This feature plays a very important role in esti- 
mating the amount of injection energy of protons for explaining the observed gamma-rays 
and the annihilation emission. 

As we see from our analysis several characteristic times determines the conditions under 
which the annihilation flux is formed. These times are: 

• The current time from the moment of eruption t: 
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• The time for nuclear p — p reactions (the hfetime of protons) Tp — nVpUpp. For gas 
densities ranging from 1 to 100 cm~^ it ranges from lO^^s to lO^^s; 

• The coohng time of positrons Tc(the thermahzation time) 

[ (IE 

where {dE/dt)c is the energy coohng rate, which will be discussed in section 7, E'e+and 
E-th are the injected energy and the thermahzation energy respectively. For a density 
n we estimate that Tc ~ 3 x 10^^(n/l cm~^)~^ s. 



In the following sections, we investigate models which are characterized by different 
ratios between these time scales. 



5. Gamma Rays from 7r°-Decay 

Pions are produced m p — p collisions whose cross-section is well-known (see, e.g., 
Stephens and Badhwar 1981; Atoyan 1992; Dermer 1986; Dogiel and Sharov 1990; Strong 
et al. 2000). 

The production function of gamma-rays is 

q^o = 7] / N{Ep)vpnHda{E, Ep) (17) 

J Ep 

where is a factor depending on the chemical composition of cosmic rays and of the inter- 
stellar medium, which for the Galaxy equals 1.6 (Stephens and Badhwar 1981), Vp is the 
velocity of protons and da is a differential cross-section for the production of a 7r° of energy 
E by a cosmic ray proton. 

According to Fatuzzo and Melia (2003) the cross-section is 

da{Ep,E^o) ctq 



dE^^o E^o 



fAx) (18) 



where x — E^o/Ep, (Jq = 32 mbarn and fT^o{x) — 0.67(1 — x)^-^ -\- 0.5e It is important 
to note that the above formula is based on the isobar model of Dermer (1986), which is 
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technically appropriate for a modest range of energies above pion decay threshold. At the 
highest energies detected by EGRET, the scahng model is more appropriate and this is also 
given in Dermer (1986). 

The emissivity of the photons produced by 7r°-decay (q^) can be calculated from (Stephens 
and Badhwar (1981)): 

oo 

q,{E„t) = 2 [ -^0^^dE, (19) 



TT 



where and E-,^ are the energy of the emitted photon and the decaying pion respectively, 

o A 

E ■ = E 



The spectrum of expected gamma-ray emission is shown in Fig. 3 together with the 
EGRET data. From the observational data we concluded that at the present time the 
constant Aq is 

= 4 X lO^^GeV^o-^ (20) 
and the total energy of accelerated protons emitting gamma-rays is 

oo 

Wp^ J EpNp{Ep)dEp ~ 1.5 X lO^Vg , (21) 

Ep 

for a density n — 1 cm~^. If we take the average density of the gas in the region of gamma-ray 
emission n — 10 cm~^ we obtain W ~ 10^^ erg for the same production rate of gamma-rays. 

Below we determine whether these protons are able to generate enough positrons in 
order to produce the observed flux of the annihilation line 



6. Production Function of Secondeiry Electrons 

Inelastic p—p collisions produce two charged pions for every neutral pion. These charged 
pions quickly decay into muons, which in turn decay into positrons and electrons, with a 
resulting emissivity 

Vmr~ " 



The integration limits in the above expressions are determined through kinematic consider- 
ations: for example for two-body decay the minimum and maximum energy in units of the 
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rest mass of the created particles are 7 — p and 7 + p, where 7 is the Lorenz-f actor of the 
source particle and p = ^^7^ — 1 is its dimcnsionlcss momentum. The lepton distribution 
from a decaying muon is given by the three-body decay probability 

dP_ ^ 8pc f u{uX-myy/' ( 47,^ \ / _ E,{E,-u) \ 

dE, P^mlc^J {pc-E, + uy \ m^cVV P'c' J ^ ' 

where u = {Ee — P^pc cos 9), p is the electron momentum and 7 is the Lorenz factor for the 
designated particle (Fatuzzo and Melia 2003). 

The production function of vr^ secondary electrons is shown in Fig. 4. The secondary vr^ 
electrons are generated above the threshold energy E >30 MeV. Below this energy secondary 
electrons are produced by Coulomb collisions (knock-on electrons) . The cross-section of KO 
electron production is (see, e.g., Hayakawa 1964; Ginzburg 1989) 

a,,c^2n{^) 24 
\mc'^ J Ej 

The composition of the production flux is an essential function of particle energy. The 
electron flux at energies above ~ 30 MeV consists of almost equal parts of positrons and 
electrons. Below this energy the contribution of positrons is negligible. 

From Eqs. (12), (20) and (22) we can estimate the total production rate of relativistic 
positrons at GC at the present moment directly associated with the current gamma-ray 
production. It equals 

'''^"-lO^V^ (25) 



dt 

Below we estimate the number of thermal positrons supplied via p — p collisions after 
their cooling by Coulomb interactions (ionization losses). The total electron production 
function for tt"*^ and KO electrons is shown in Fig. 4 where we used the results of Daniel and 
Stephens (1975). 

The distribution function of these relativistic positrons is described by Eq.(9). In the 
relativistic energy range electrons lose energy by synchrotron and inverse Compton emission 

' ^ ' Wph + Q.2lO''—)E^ = -/3iE^ (26) 



dt 9 \mc^ J mc^ \ Stt / 

The characteristic time of energy loss r in this case is 



E 
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Then for the proton spectrum (10) and the production function (22) the spectrum of sec- 
ondary electrons can be derived from the general equation (9), see Syrovatskii (1959) and 
Gratton (1972), 



I oo 



exp 



(28) 



X exp 



[r - ro 



\2-\ 



S{t-to-T) 



where 



(29) 



Qe is derived from Eq. (22), and 9{x) is the Heaviside function. 



If there were several eruption processes we should sum over them in order to get a 
resulting spectrum of positrons. 



7. Evolution of the Positron Spectrum 

7.1. Kinetic Equation 

In order to estimate the flux in the e"*" — e~ annihilation line we should calculate how 
many thermal positrons can be supplied by p — p nuclear reactions under the influence of 
energy losses. Similar calculations were performed by Stecker (1969) who analyzed effects 
of energy losses only. Crannell et al. (1976) calculated the positron density produced in 
solar flares from the Fokker-Planck equation when positrons are injected at an energy of 
1 MeV. Recently Guessoum et al. (2005) have studied the lives and deaths of positrons 
in the interstellar medium by using Monte Carlo simulations as well as a binary reaction 
rate approach. In our case we should use the generalized Fokker-Planck equation because 
positrons are injected at a relativistic energy of about 30 MeV. 

The generalized Fokker-Planck equation, describing the particle distribution function / 
written in dimensionless variables is (see Dogiel 2000; Liang et al. 2002) 



dt uq 



(30) 



where Qg describes the distribution of sources emitting fast particles. If we compare this 
equations with Eq.(9) for relativistic and subrelativistic cosmic rays we see that in order 
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to describe evolution of the positron spectrum from relativistic to thermal energies we in- 
clude the terms in rhs which are responsible for the formation by Coulomb collisions of the 
equilibrium Maxwellian distribution with a given temperature of the background gas. The 
injection spectrum is given by the term q^. which is determined by the production of sec- 
ondary electrons at energies above 30 MeV. As one can see from Fig. 5 the change in / 
at thermal energies implied by Eq. (30) is then a simple scaling with time, mediated by 
the ionization loss timescale. The charge-exchange annihilation processes describing by the 
cross-section act for positrons as an absorbing boundary at low energies, especially for 
eV temperatures when the cross-section of annihilation increases by several orders of mag- 
nitude (Bussard et al. 1979; Guessoum et al. 2005). Therefore, for a constant production of 
positrons we have in the stationary case the annihilation rate of thermal positrons exactly 
equal to J qe{p)p^dp. 

In our case this source function describes the production of positrons (shown in Fig. 
4) with the rate ~ 10"^^ e"*" s~-'^GeV~-'^, the main part of which is concentrated near energies 
Ee+ ~ 30 — 50 MeV, and processes of positron annihilation are described by cross-sections 
(Tee s-nd (7 an (sec sectiou 6.1). 

We define non-dimensional variables p — pj \JmkT and t — i-UQ, where the characteristic 
frequency vo is 

2'KneC^rlme 

J^o - — , , i3ij 

where Ue is the average electron number density, and is the classical electron radius 



2 / 2 
e /TTleC 



lcm-3 and = 2.5 eV, z/q ^ 6.7 x lO^^^ s^^ 



In Eq. (30), B describes the particle energy gain and loss and A the momentum diffusion 
due to Coulomb collisions with background particles. 



A{p) = p" 



dp 
It 



7 



(32) 



B{p) = p' 



dp 



dp 



dt 



synIC 



where ^ — 1 + E{p)/meC^ is the Lorentz factor. 



dp 
'di 



brem 



(33) 



The dp/dt terms are rates of momentum loss due to ionization, synchrotron plus in- 
verse Compton (IC) and bremsstrahlung emission respectively (e.g. Hayakawa 1969). The 
ionization loss term is 



dp 
'di 



'pY ^ kT, a/t^ 



X 



log 



£(j;)/;;,r^(-^-l) 



+ 0.43 



(34) 
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where the plasma frequency Up — y^Ane'^ne/me- The combined synchrotron and IC loss term 
is 



/ synIC 



321 



rrieC^ 1 



X 



nrt 



= [Ucmb + Umagl I 



E{p) 



(35) 



where Ucmb is the energy density of background photons, and Umag is the magnetic energy 
density. The bremsstrahlung loss term is given by 



dp 
'dt 



brem 



meC^ 1 E(p) ngmgC^e^r^ 
kT^ vq rUeC^ hc\/mekTx 



log 



2{E{p) ^ rrieC^) 1 



(36) 



The corresponding times of synchrotron+inverse Compton energy losses are shown in Fig. 
1. The magnetic field strength is if = 5 x lO^^G and the energy density of background 
photons Ucmb = 1 eV cm"^. 



7.2. Numerical Calculations of Thermal Positron Numbers for the Case of 

Quasi-Stationary Production {t <Tp ^ Tc) 

For calculations of the positron number at different times we used Eq. (30) describing 
the evolution of the positron spectrum in both the relativistic and nonrelativistic energy 
ranges. The dimensionless time t — 1 corresponds to 1.5 x 10^^ s. We include into the 
equation the process of positron annihilation whose cross-section at different energies is (see, 
e.g., Crannell et al. 1976) 

nrl [f + Aj+l r r^-rl 7 + 3 \ 
= 7TI [ 7^-1 b+ VT^J - -^^^ j (37) 

If the energy of positrons is below 1 keV and the medium is weakly ionized (as is the case 
at temperatures below 10 eV) then we should also include the charge-exchange annihilation, 
whose cross-section is : 



<7ee = 1.2 X 10-^^ • - - ^) 

where Ef,+ is the positron energy and Eq — 6.8 eV. This equation was obtained from the 
experimental data of Sperber et al. (1992) with the approximation of Fatuzzo et al. (2001). 
The corresponding time of positron annihilation at different energies Tan — i^cTan)'^ is 
shown in Fig. 1 by a dotted line. 
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The production function of positrons is shown in Fig. 4. The injection function can 

be presented as qe{p,t) oc Qo{t)9{p — Po)p~^ where po corresponds to the maximum of the 
function qe{p) at the energy E = 30 MeV. We assume that annihilation processes take 
place in the extended medium surrounding the region, where the injected protons produce 
relativistic positrons. The time variation of the injection function is caused by catastrophic 
and ionization losses of protons. However, if we consider a time t < Tp we can neglect time 
variations of the production function corresponding to the almost stationary case of positron 
production. We normalized the distribution function 

oo 

J P^dp q,{p) = Qo (39) 



where Qq can be obtained from Fig. 4 and is about 10"^^ s~^ for t < Tp. 

If we calculate from the kinetic equation (30) the number of thermal positrons anni- 
hilating in the cV temperature medium on the assumption that they are produced by the 
same relativistic protons which are responsible for the observed gamma-ray flux of 2 ■ 10^^ 
erg s~^ then we obtain the value 2 ■ 10^^e"'"s~^, which is two orders of magnitude less than 
necessary to explain the annihilation flux from the Galactic center. Therefore the stationary 
model fails to explain the data. Thus, in the framework of this model we cannot satisfy both 
conditions 1) and 5) presented at the end of Section 2. 



8. Nonstationciry Model of Positron Production 

8.1. Positron Spectrum {t > Tp ^ Tc) 

One of the solutions to this problem might be that in the past the positron production 
rate was much higher than follows from the current gamma-ray flux from the GC. Under 
this condition the production function of positrons calculated from Eq.(12) is an essential 
function of time, q{p,t). 

In Fig. 5 the evolution of the positron spectrum calculated from Eq. (30) is shown for the 
parameters of the interstellar medium: the density n — 1 cm~^, the temperature T — 2.5 eV, 
the magnetic field strength H — 5 x 10~^G and the energy density of background photons 
Ucmb = 1 eV cm"2. 

The distribution function of the positrons injected at the initial stage at p > po is 
shifted to the momentum region p < po under the influence of Coulomb collisions that can 
be seen in Fig. 5 as a "bunch" of positrons propagating into the region of small momenta. 
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Later Coulomb collisions start to form the equilibrium Maxwellian distribution in the range 
E < 10~^ keV . At a subsequent point of time Coulomb collisions continue to form the 
equilibrium distribution and accumulate positrons in the thermal energy range, as can be 
seen in Fig. 5 as the increasing Maxwellian distribution. At the final stage when the source 
ceases to work the number of thermal positrons decreases because of annihilation (see Fig. 
6). 

We notice that the positron evolution and the number of thermal positrons depend on 
the background gas temperature, e.g. in the plasma with the temperature 100 cV we can 
produce a much larger amount of positrons because the annihilation cross-section in this 
medium is several orders of magnitude smaller than in a neutral gas (see Fig. 7). 

Now we proceed to calculate the intensity of emission due to annihilation of positrons. 
Depending on the temperature and density of the ambient medium, electrons and positrons 
can annihilate either directly or via formation of positronium (Bussard et al. 1979; Guessoum 
et al. 1991, 2005). 

The in-flight differential spectrum of the 7-rays produced by annihilation of a positron 
on the ambient electrons with density Ug (Aharonian and Atoyan 1981; Aharonyan and 
Atoyan 2000): 

qan{e) = — + — + (40) 



7+P+ 



7+ + 1 - £ e 



1 1 \ /I 1 



e 7+ + 1 - \E 7+ + 1 - £ 



where 7+ = E^/nieC^ is the Lorentz-factor of the positron, = a/7+ — 1 is the dimension- 
less momentum of positron and e = E /rrieC^ is the dimensionless photon energy. The energy 
spectrum of in-flight annihilation 



Nan{,e) = / qan{no,e)ne{e)d^Q (41) 



7 



The lower integration limit can be obtained from the kinetic equations and equals 

,ie) ^ (42) 

The energy spectrum of annihilation photons produced by the charge-exchange process 
is (Bussard et al. 1979) 

N,,{e) = J dEonHne{Eo)ace{Eo)v - [{Eo - 6.SeV){Eo - 6.SeV + Ae)]'/^ (43) 



E 
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where E — [e — SqY /e-\- 6.8ey, Sq - the rest mass of positronium, uh - the density of neutral 
hydrogen. 

To fit the INTEGRAL data (Churazov et al. 2005), shown at Fig. 8, we used the 
following medium parameters: plasma temperature 2.5 eV, plasma density 1 cm~^, degree 
of ionization 5% (see, e.g., Kaplan and Pikelner 1979). The plasma temperature affects the 
width of the annihilation line. In order to fit the INTEGRAL data for energies at and below 
511 keV we need a temperature in the region 2-5 eV. 

The line emission at these temperatures is due to the two photon decay of the charge- 
exchange process. The emission below 511 keV is generated by three photon decay of positro- 
nium. The process of in-fiight annihilation in a low temperature medium is negligible be- 
cause the time for thermalization in this medium is smaller than the characteristic time of 
the in-fiight annihilation. This process is essential in the hot component of the interstellar 
medium, about 20% of which is filled by a hot plasma with the parameters T = 100 eV and 
n = 3 X 10~^ cm~^. However, as one can see from Fig. 9 the flux of the annihilation line 
produced in the hot interstellar plasma by in-flight annihilation is negligible in comparison 
with emission in a cold plasma. 

The temporal variations of the gamma-ray and the annihilation fluxes are shown in Fig. 
10. The maximum of the fluxes are taken as unity. We see that the gamma-ray flux is 
always decreasing with time, while the annihilation flux reaches its maximum a long time 
after the proton eruption (~ 3 x 10^^ s) and then gradually decreases. We can estimate 
the necessary input proton energy without going into the technical calculation details. The 
current observed energy flux of gamma-rays is 2 x lO^^erg s~^ and the current annihilation 
rate of thermal positrons is 10"^"^ positron s^^, corresponds to a higher proton collision rate 
in the past. From Fig. 10, one can see that this situation is realized if we are observing 
these fluxes at a time t > Tp when the current gamma-ray intensity is at least three order of 
magnitude less than its maximum, while the current annihilation flux falls from its maximum 
value by less than one order of magnitude (just because of the delay effect). In this case, 
the necessary input proton energy is at least ~ lO^^erg. Of course, this estimate sensitively 
depends on the number density. Even if we take the average density of the gas in the bulge 
as high as n = 100 cm~^ (see Section 4) we get a value of the energy necessary to produce 
the gamma-ray and annihilation fluxes about Wp ~ 10^^ erg. This amount of energy may 
be able to be supplied by a very massive star accreting onto the black hole (see Eq. 5) or by 
extracting the rotation energy of the black hole (see Eq. 6). However, we may not be able 
to obtain a distribution of the annihilation line extending to several hundred pc from the 
GC because the characteristic length of diffusion propagation is about 100 pc for this high 
gas density. Thus, in the framework of this model we have problems with the conditions 2) 
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or 4) from Section 2. Nevertheless, we conclude that this model is marginally acceptable if 
a 50 Mq massive star was captured by the black hole about hundred million years ago. 

9. Spatially Non-Uniform Model (r^ < i ~ Tc) 

To solve the possible problem of energy excess we take into account that the background 
gas is heterogeneously distributed in the bulge (see Section 4). According to Jean et al. 
(2005) hydrogen gas in the nuclear bulge is trapped in small high density molecular clouds. 
Penetration of cosmic rays into molecular clouds was analyzed in a number of papers (see, 
e.g., Cesarsky and Voelk 1978; Dogiel and Sharov 1985; Padoan and Scalo 2005). From 
the EGRET observations it is known that GeV cosmic ray protons freely penetrate into 
molecular clouds (Digel et al. 1999) while it is almost unknown how MeV electrons interact 
with the clouds. However, Skilling and Strong (1976) showed that the flux of cosmic rays 
below a few hundred MeV might be completely excluded from molecular clouds. 

In our case it means that relativistic protons propagate in the medium with the average 
gas density about n=30 - 1000 cm~^ (depending on their propagation distance from the bulge 
center) while the secondary positrons ejected from the clouds propagate in the intercloud 
medium only, where the average gas density is about n ~ 1 — 10 cm"^. 

For protons it means that they fill the region of a radius < 100 pc (if the diffusion 
coefficient equals the average in the Disk, D ~ 10^^ cm^s~^) during their lifetime Tp ~ 
(1 — 5) X 10^^ s while electrons can fill a sphere of radius ~ 400 pc over the period of time 
of their Coulomb cooling (~ 3 x 10^^ s for n = 1 cm""^). 

The time variations of the gamma-ray and annihilation emission are shown in Fig. 1 1 for 
an average density for protons n = 30 cm^^. As is to be expected the gamma-ray flux drops 
rapidly from its initial value while the annihilation flux reaches its maximum ~ 3 x 10^^ s 
after the eruption. The energy of relativistic protons necessary to get the peak production of 
thermal positrons of about 10^^ e+s~^ is about 7-10^^ erg. The energy in relativistic protons 
changes weakly if we increase further the gas density. Thus, for a density n — 100 cm~^ the 
necessary energy of the protons is 5 ■ 10^^ erg. These injected proton energies are compatible 
with the estimates in Eq.(5) or Eq.(6) for a massive star capture. 

Besides, if the average density of the gas in the medium traversed by positrons is a 
little bit higher than 1 cm~^ (e.g, 3 cm~^) than we obtain slightly smaller values for the 
thermolization time (~ lO^^s) as well as the necessary energy output of protons (~ 10^^ erg) 
that makes our estimates even more acceptable. 
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Then, we can assume that the gamma-ray and the annihilation fluxes from the Galactic 
center belong to different accretion events. Gamma-rays were generated relatively recently 
(~ 10^^ s ago) when a one solar mass star was captured by the black hole. Indeed the 
necessary amount of energy in injected protons is about > 10^^ — 10^^ erg (see Section 5). In 
this respect it is clear why this emission was observed by EGRET as point-like. The average 
length of proton propagation is < 100 pc for the diffusion coefficient D = 10^^ cm^s~^. If 
the annihilation line emission was produced ~ 10^^ s ago by a previous event of a 30 Mq 
star capture, with an energy release of > 10^^ erg in protons, then it appears as an extended 
source long after the time when the gamma-ray emission produced by this capture died out. 
Thus, this model satisfies the conditions formulated in Section 2 though the annihilation and 
gamma-ray fluxes observed from the GG are not connected with each other in this case. 



In the framework of the non-uniform model we can estimate the minimum energy at 
relativistic protons in order to produce in the current time the thermal positron rate of about 
10^^e"'"s~^ and not contradict the gamma-ray observations. We remember that the present 
gamma-ray flux from GC is about 2 ■ 10^^ erg ~^ and it is seen as a point-like object. We 
remember also that the gas distribution is unknown but that the main part of the gas is 
concentrated inside the central region of the bulge where the its density is as high as ~ 1000 
cm~^ (see Jean et al. 2005). If relativistic protons freely penetrate into the clouds the 
average density n ~ 1000 cm~ . We define this density 

The diffuse gas is distributed in the intercloud medium with an average density of 10 
cm~^ and it may reach at the central regions the value of 30 — 60 cm^'^. If positrons do not 
penetrate into the cloud then the gas density of the medium in which positrons are cooled 
down by the Goulomb collisions is determined by these values. Below, we define this gas 
density as rice- Let us estimate the minimum energy of primary protons necessary to produce 
the annihilation fiux as a function of Upp and rice- 

In the framework of the nonuniform model the initial energy of protons necessary to 
produce the annihilation emission depending on Upp and rice is varying in the following way: 
It is clearly seen that when we increase the value of Upp the output energy of protons decreases 
because the production rate of relativistic positrons is proportional to the gas density 



However as follows from Table 1 the proportionality between Wp and Upp is not linear as 
it would seem from Eq.(44). Prom Eq.(44) one can see that the higher the density Upp the 



10. Continuous stellar captures 




(44) 



-23- 



Table 1: Injected Proton Energy 



npp{cm 3) 



ncc(cm 2) 
1 



W (erg) 
^055 

7 • 10^3 
5 ■ 10^3 
2 ■ 10^3 

2 ■ 10^3 

1.5 ■ 10^^ 
4.5 • 10^2 

3 • 10^2 



1 

30 

100 
100 
1000 
1000 
1000 
1000 



10 

1 

10 
30 
60 



1 



1 



shorter is the injection impulse of positrons which is determined by the hfetime of protons 
Tp = nppGppC. The necessary energy output of protons as a function of the ratio Upp/ricc is 
shown in Fig. 12. In the figure we take the density Ucc — 1 cm~^ and Wq is the proton energy 
when Upp/ricc — 1- The variation of the proton injected energy Wp can easily be understood 
from the solution of the diffusion equation. When Upp/n,.,. = 1 we have Tp > Tc and the 
energy Wp at this density for a fixed flux of thermal positrons at the level Fe+ ~ lO'^'^e^s""^ 
is of the order of Wp oc Tp ■ i.e it decreases as 1^ oc l/upp. However when Upp ^ 
the impulse of proton injection is very short Tp ^ Tc (almost delta-function injection) and 
we obtain Wp oc Tc ■ -F+, i.e Wp is independent of Upp as one clearly sees in Fig. 12. The 
time Tc cannot be too short otherwise the gamma-ray flux would not drop down below the 
observed gamma-ray flux from the GC (see Fig. 11). Variations of the gamma-ray flux for a 
flxed annihilation flux of thermal positrons of 10^^ ph s"^ as a function of the ratio Upp/ricc 
is shown in Fig. 13. From this flgure we see that this ratio cannot be below Upp/ricc = 17 
otherwise the peaks of gamma-rays and the annihilation emission are so close to each other 
that the flux of gamma-rays does not have enough time to drop below the observed limit. 

From this restriction we conclude that the injected proton energy cannot be less than 
~ 3 X 10^^ erg even for Upp = 1000 cm^"^ while Ucc = 60 cm~^. This gives the average distance 
of proton propagation by difi^usion of about 20 pc and about 100 pc for positrons during their 
lifetime, which is less than necessary, but the diffusion coefficient may be spatially dependent 
such that it gives a more extended positron spatial distribution. Under this assumption even 
the capture of 1 Mq star can give the necessary energy in protons as follows from Eq.(5) or 
Eq.(6), and the process of annihilation emission is not an exceptional event as we had in the 
previous models when a massive star capture was needed to supply the necessary energy in 
relativistic protons. 
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Since there are so many stars at the GC, there is very good reason to beheve that the 
capture of stars can continuously take place from time to time. If there are several successive 
eruptions near the black hole then the annihilation flux is fluctuating near a certain level of 
emission as is shown in Fig. 14. In the calculations we take n^c = 60 cm~^ and Upp = 1000 
cm~^ which give Tc ~ 5 x 10^^ s and ~ 7 x 10^^ s respectively. The time between 
two successive eruption equals ~ 8 x 10^^ s. When the variations of the energy fluxes of 
gamma-rays and the photon flux of the annihilation emission are normalized relative to their 
maximum their values, which for gamma-rays is 3 x 10^^ erg s~^ and for the annihilation 
flux is 10^^ ph s~^, the average injected proton energy is ~ 8 x 10^^ erg. 



11. TeV Gamma-Rays 

It is very interesting to find other confirmations of eruption processes from observed 
data. In this respect, the new data in the TeV gamma-ray range are of strong interest. 

Recently, TcV 7-ray emission at a luminosity higher than 10^^ ergs s~^ from the direction 
of the galactic center has been reported by three independent groups, Whipple (Kosack et al. 
2004), CANGAROO (Tsuchiya et al. 2004), and HESS (Aharonian et al. 2004). The most 
plausible candidates suggested for this emission include the black hole Sgr A* (Aharonian 
and Neronov 2005b) and the compact and powerful young supernova remnant (SNR) Sgr A 
East (Crocker et al. 2005). The angular scale of the TeV source was determined by HESS 
to be about a few arc- minutes, indicating that this 7-ray source is located in the central 
~ 10 pc region (Aharonian et al. 2004). The gamma-ray spectrum observed by HESS can be 
described by F{E^) = (2.5 ± 0.21) x 10~^^i?~^-^^ ph/{cm'^ sTeV), which is consistent with 
decay of 7r° produced by p-p collisions with the initial spectrum of the injection protons 
as Q{E) = QqE'"^-^ exp{—E / Eq) , where the cutoff energy is Eq = 10^^ eV (Aharonian and 
Neronov 2005b). In comparing with the diffuse gamma-ray spectrum detected by EGRET 
(cf. Eq. 2), the TeV spectrum is flatter. If both GeV and TeV gamma-rays result from the 
decay of 7r° which are produced by p-p collisions, it appears that these two different energy 
gamma-ray spectra must be produced by two different groups of protons. Furthermore, 
the protons responsible for the diffuse GeV gamma-ray emission must be ejected by Sgr A* 
much earlier than those responsible for the TeV gamma-rays because their angular sizes are 
so much different. In fact, these two requirements are consistent with our model considered 
here, i.e. the galactic black hole is activated by capturing stars. 

If we assume that positrons arc produced via p — p collisions and these protons arc emit- 
ted from the galactic black hole, by taking the typical diffusion coefficient (D~ lO^^cm^s^^), 
we can estimate that these relativistic protons should be emitted at a earlier time t^iff ~ 
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6xio^W^/s ~ 3 X 10^ s, which is consistent with the requirement of the annihilation hne. 
Again if we assume that TeV photons also result from the decay of neutral pions produced 
by p-p collisions and that these protons are emitted from the galactic black hole, then taking 
the same diffusion coefficient, we obtain that these protons are emitted at a earlier time 
idiff ~ 6xio°Vm2/s ~ 5 X lO^ycars, which is shorter than the typical capture time of a main 
sequence star by a factor of 60. However, D~ 10'^^cm^s~^ is the mean diffusion coefficient 
for the Galaxy obtained by cosmic ray measurements (Berezinskii et al. 1990). The dif- 
fusion coefficient near the Galactic center can be significantly smaller. For example, the 
Bohm diffusion time can be estimated as tBohm ~ (^)^ ic '^^^^re A is the diffusion distance, 
Ep/eB is the Larmor radius, Ep is the proton energy for producing TeV photons and 
B is the local magnetic field. Using A ~10 pc, Ep ~3TcV and B ~ 10~^G, we obtain 
tBohm ~ 2 X lO^years, which is similar to the requirement of the continuous capture model. 
Taking the observed TeV luminosity ~ lO^^erg/s that implies the injected proton energy 
must be at least ~ lO^^crg (Aharonian ct al. 2004), which is consistent with our estimation 
of energy release for capturing a IMq star (cf. Eq. 5). Furthermore, the injected proton 
index for TeV gamma-rays is ~ 2.2 but the injected proton index for the diffuse GeV gamma- 
rays is ~ 2.6 in our model calculation. It appears that they have different origins. However, 
this discrepancy can be explained by the following possibilities. (1) The simplest explanation 
may be that the injected proton spectrum can vary from injection to injection. (2) The most 
realistic explanation should be due to propagation effects and leakage effects which occur in 
the Galaxy where the injection spectrum with the spectral index as Q{E) oc _E'~(2.o-2.2) g|;ggp_ 
ens due to energy variations of the spatial diffusion coefficient, as D oc E~^'^ (see Berezinskii 
et al. 1990). 

Lu et al. (2005) have suggested that the TeV emission from the GC resulting from a 
powerful relativistic proton jet emitted by Sgr A* about ~ lO'^ycars ago, which result from 
the capture of a red giant star. However, their conclusion is based on a comparison of the 
diffusion length of the relativistic protons with the upper limit on the size of the TeV source. 
It should be noticed that the lack of information beyond 10 pc does not mean an absence of 
such an extended component. Detection of an extended region with a weaker signal is very 
difficult. In this model we predict that when more TeV data are available the TeV spectrum 
may show a spatial dependence. In particular, the TeV signals from a more extended region 
should exhibit a steeper spectrum if the propagation and leakage effects are indeed the origins 
of the spectral discrepancy between GeV gamma-rays and TeV gamma-rays. 
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12. Discussion 

We suggest that the capture of stars by the Galactic black hole is a natural energy 
mechanism for providing relativistic protons with a typical energy of ~ lO^^^^^erg. An 
essential point of our analysis is the non-stationary processes and the non-uniform density 
conditions of emission regions. 

The duration of energy emission as well as of the existence of relativistic protons af- 
ter accretion is much shorter than the mean interval between two accretion events. Once 
ejected by the central black hole, relativistic protons propagate by diffusion into the inter- 
stellar medium surrounding the hole occupying a more and more extended region around 
the Galactic Genter. GoUisions of these protons with the ambient gas result in secondary 
gamma-rays and relativistic positrons with energies £■ > 30 MeV. The total gamma-ray 
flux is almost constant during the lifetime of the protons but then it decays exponentionally 
with the characteristic time dependent on the gas density. Secondary relativistic positrons 
produced immediately via p — p collisions are almost unseen at the first stage after the erup- 
tion because the annihilation fiux of fast positrons is almost negligible. These positrons are 
cooled down by Goulomb collisions with the ambient gas with the characteristic time for 
Goulomb collisions. Only when positrons reach an energy of about several eV they effec- 
tively annihilate in the warm neutral and ionized phases of the interstellar medium where 
the cross-section for annihilation increases by several orders of magnitude (see Bussard et 
al. 1979; Guessoum et al. 1991, 2005) and results in an increase of the annihilation emis- 
sion. The maximum of the annihilation emission is reached long after the eruption moment. 
At that time positrons fill an extended region around the Galactic center by diffusion and 
therefore in all cases the annihilation emission can be seen as an extended region around 
the Galactic center, the annihilation flux is never seen as a point-like central source. On 
the other hand, the high energy gamma-rays are produced from the very beginning of the 
eruption and can be seen in the first stages after the eruption as a point-like source. 

Our calculations show that it is marginally possible that the gamma-rays coming from 
the GC and the annihilation emission belong to the same eruption process provided the 
injected proton energy is ~ 10^^ erg, this might be possible theoretically if we take model 
parameters to their extreme values. An alternative assumption that the currently observed 
gamma-rays from the GG and the currently observed annihilation emission belong to at least 
two different injections seems to be more attractive. The first one belong to a recent capture 
of a one solar mass star with the injected energy of about 3 x 10^^ erg while the second one 
is a consequence of a previous capture of a much more massive star with an energy release as 
high as 10^^ erg. This does not contradict estimations from theoretical models if we assume 
capture of a massive star. 
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Purthermore, if processes of secondary particle production occur very close to the GC 
where the gas density is about 1000 cm~^ then even a 1 Mq star capture can provide the 
necessary energy of protons, which is in this case ~ 3 ■ 10^^ erg. The GC is a region with 
a high concentration of stars, the capture of a star should continue to take place from time 
to time. Although the characteristic capture time scale is not known, TeV emission from 
the GC suggests that it should be > lO^yr. If this is true, assuming the positrons cannot 
penetrate the cloud, so that the cooling time is longer than the p — p collision time, our 
calculations indicate that the positron annihilation rate is more or less constant whereas the 
current gamma-ray flux is measured somewhere between two proton injection events. This 
naturally explains why the annihilation rate is about two orders of magnitude higher than 
follows from the current emission rate of gamma-ray photons. 
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Fig. 1. — Variations of life-times due to ionization losses when n — 1 cm"^: for electrons 
(dashed-dotted line) and protons (solid line), due to synchrotron {H — 5 x 10"^ G) and 
inverse Compton {wph = 1 eV cm~^) energy losses for relativistic electrons (dashed line), 
and due to direct annihilation (dotted line). 
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Fig. 2. — The spectrum of protons necessary at different moments of time for an ambient 
plasma density n — 1 cm~^. 
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Fig. 3. — The flux of gamma-ray emission in the direction of the Galactic center measured 
with EGRET and the results of numerical calculations. 



-36- 



10 




-| Q'^'^ I , , I , , I , , I , , I , , I 

— -9 — 1 n 1 9 

10 10 10 10 10 10 

E (GeV) 



Fig. 4. — The total production function of secondary electrons from 7r^-decay and Coulomb 
collisions (solid line) and the production function of positrons (dotted line). 
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Fig. 5. — Time variation of the positron distribution function f{p) for a gas temperature 
T = 2.5 eV. The dimensionless momentum is defined as p/\^rn^kT so the dimensionless 
momentum at unity corresponds to a positron with energy ~ kT. 
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Fig. 7. — Time variation of the amount of thermal positrons produced by relativistic protons 
injected 4 x 10^^ s ago in the gas with a temperature T — 100 eV and a density n = 3 x 10~^ 
cm~^. 




Fig. 8. — The total annihilation spectrum from a plasma with a temperatures 2.5 eV, together 
with the INTEGRAL data 
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Fig. 9. — The annihilation flux from the hot plasma of the central region of radius ~ 0.7 kpc 
with a temperature of 100 eV, a density of 3 x 10""^ cm~^ and a flUing factor of 30%. 
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Fig. 10. — Time variation of the gamma-ray (dashed hne) and annihilation (sohd hne) fluxes 
for a gas density n — 1 cm~^ 
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Fig. 14. — Gamma-ray and annihilation emission from several successive eruption of protons 
by the black hole. The parameters are presented in the text 



